Non-Gaussianity and the CMB Bispectrum: confusion between primordial and 
primary-Lensing-Rees Sciama contribution? 
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We revisit the predictions for the expected Cosmic Microwave Background bispectrum signal 
from the cross-correlation of the primary- lensing-Rees-Sciama signal; we point out that it can be a 
significant contaminant to the bispectrum signal from primordial non-Gaussianity of the local type. 
This non-Gaussianity, usually parameterized by the non-Gaussian parameter Jnl, arises for example 
in multi-field inflation. In particular both signals are frequency- independent, and are maximized for 
nearly squeezed configurations. While their detailed scale-dependence and harmonic imprints are 
different for generic bispectrum shapes, we show that, if not included in the modeling, the primary- 
lensing-Rees-Sciama contribution yields an effective /nl of 10 when using a bispectrum estimator 
optimized for local non-Gaussianity. Considering that expected 1-cr errors on /ml are < 10 from 
forthcoming experiments, we conclude that the contribution from this signal must be included in 
future constraints on }nl from the Cosmic Microwave Background bispectrum. 

PACS numbers: 



I. INTRODUCTION 

The increased sensitivity of the forthcoming Cosmic 
Microwave Background (CMB) experiments will open 
the possibility to detect higher-order correlations in the 
CMB temperature fluctuations beyond the power spec- 
trum. This means that it would be possible to study 
in detail eventual deviations from Gaussian initial condi- 
tions and thus gain an unique insight into the physics 
of the early universe (see e.g.[l| for a review). Since 
gravitationally-induced non-linearities at last scattering 
are much smaller than in the late-time universe, the CMB 
is expected to be the best probe of the primordial non- 
Gaussianity (e.g., @,S]). 

Moreover, even in absence of these primordial devia- 
tions, measuring the CMB three-point correlation func- 
tion or, equivalently its Fourier analogue, the angu- 
lar bispectrum, would be very useful to trace the im- 
print of the non-linear growth of structures on secondary 
(i.e. late-time) anisotropies and would open a new win- 
dow into the understanding of the evolution and growth 
of structures. The expected bispectrum signature of 
secondary CMB anisotropies has been studied in e.g., 
0, H H 0, @, i, [T(J [Tl|]. In addition, non-linear physics 
happening between the end of inflation and the last scat- 
tering surface may leave some imprints in the CMB bis- 
pectrum (e.g., [2, [H, [H, EH E3 for local- type non- 
Gaussianity, which is relevant here). 

To clarify the use of our nomenclature, by primary 
non-Gaussianity (or primary bispectrum) we refer to the 
combined effect of primordial non-gaussianity and of the 
physics happening between the end of inflation and the 
last scattering surface. This name is chosen accordingly 
to the CMB nomenclature of "primary anisotropies" 
which are related to the primordial ones but are further 
processed. In the literature what we call "primary" non- 
Gaussianity is often loosely referred to as "primordial" . 



Analogously, we use secondary non-gaussianity (or sec- 
ondary bispectrum) to refer to late-time physics, as it 
is usually done for secondaries CMB anisotropies. Note 
that the Integrated Sachs Wolfe effect is thus considered 
a secondary anisotropy. 

After galactic foregrounds, point sources and the 
Sunyaev-Zeldovich signature from clusters are expected 
to be the dominant source of non-Gaussianity in CMB 
data. At I < 1500 the statistical properties of these two 
contributions is expected to be quite similar: both are 
expected to behave as an extra Poisson contribution (see 
[1 81 ] for treatment). In addition, both signals have a well 
known frequency dependence that can be used to clean 
CMB maps from this contaminating signal. 

In [3] it was shown that the leading contribution to the 
CMB secondary bispectrum with a blackbody frequency 
dependence is that of the primary-lensing-Rees Sciama 
correlation, where by Rees-Sciama we mean the combi- 
nation of the linear effect (Integrated Sachs Wolfe) and 
of the non-linear one. Note that in the original paper 
it was called primordial-lensing-Rees Sciama, because of 
the loose nomenclature convention explained above. Ref 
7] also computed the expected signal-to-noise for this ef- 
fect and found that experiments such a Planck 1 or ACT 
2 should be able to obtain a high statistical significance 
detection. This was then independently confirmed by 
[l]|[2(|. Note that there are sources of non gaussianity 
that are not strictly primordial (inflationary) but they 
arise between the end of inflation and the last scattering 
surface. These may well dominate over the inflationary 
contribution but they are mostly of equilateral type (e.g., 
0, EE EH ; as we will see below, here we concentrate on 



1 www.rssd.esa.int/Planck 

2 http:/ /www. physics. princeton.edu/act/ 
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the local (squeezed) type. 

Sparked by a recent study claiming a more than 95% 
confidence limit evidence 21| of local non-gaussianity, 
but see [H, US HI, H3|, the subject of primordial non- 
Gaussianity has received renewed attention. We are thus 
motivated to revisit the predictions for the expected bis- 
pectrum signal from the primary- lensing-Rees Sciama 
correlation in light of new developments since the year 
2000: a much better-determined fiducial cosmological 
model with a lower erg, improved description of linear 
and non-linear evolution of clustering in the presence of 
dark energy, optimized estimators for primordial non- 
Gaussianity via the bispectrum signal and the tantalizing 
hint of a possible detection. Instead of concentrating on 
the usefulness of the primary-lensing-Rees Sciama bis- 
pectrum in constraining dark energy as done so far in 
the literature, we will explore whether it could be con- 
fused with the primary (or primordial) non-Gaussian sig- 
nal and examine possible ways to separate the two. 

The rest of the paper is organized as follows: in Sec. |TT] 
we review the basic of the primary CMB bispectrum and 
of the primary- lensing- Rees Sciama one. In Sec. IIIII we 
present numerical results for the expected signal-to-noise, 
the dependence of the signal on bispectrum shape and 
we quantify the dependence of the secondary bispectrum 
signal on different descriptions for non-linear evolution of 
clustering. In Sec. HVI we explore a possible confusion be- 
tween the two bispectra and prospects for separating the 
signals. Finally we conclude in Sec. [V] Useful formulae 
are reported in the Appendix. 



II. THE PRIMARY AND THE SECONDARY 
LENSING-REES-SCIAMA CMB BISPECTRUM 



In this section we review the necessary background and 
the basic description of the primary CMB angular bispec- 
trum and the secondary one arising from the cross cor- 
relation among the primary-lensing-Rees Sciama (L-RS) 
contributions. For the primordial non-Gaussianity we 
will consider the so-called local type characterized by a 
momentum- independent non-Gaussian parameter -Jnl~ 
and refer to [3, la, HH . This is the workhorse model for 
testing deviations from gaussianity both in CMB data 
and in large-scale structure data. In reviewing the sec- 
ondary L-RS bispectrum we mainly follow @, • 



A. Primordial Non-Gaussanity 



bations 3 in momentum space, $(k), as: 
($(k 1 )*(k 2 )*(k 3 )> = (27r) 3 ,5 3 (k 1 + k 2 + k 3 )F(fc!, k 2 ,k 3 ), 

where the function F{k\, fc 2 , fe) describes the correla- 
tion among the modes and depends on the shape of the 
(ki, k 2 , k 3 ) triangle in momentum space. Different mod- 
els make different predictions for the function F, depend- 
ing on the mechanism of production of such a correlation 

There are two main, physically-motivated classes (32[: 

i) Local form (Squeezed configurations). 

This non-Gaussianity arises from the non-linear relation 
between the light scalar field (different from the infla- 
ton) driving the perturbations and the observed In 
the weak non-linear coupling case the non-Gaussianity 
can be parametrized in real space as in Eq. ^ with Jnl 
quantifying the "level" of non-linearity. Being local in 
position space, it couples Fourier modes far outside the 
horizon. The signal is maximal for squeezed triangle con- 
figurations with the coupling of one large-scale mode with 
two small scale modes. Examples of this class of models 
can be the curvaton scenario [33j or the Ekpyrotic model 

dm. 

ii) Non-local form (Equilateral configurations). 

In this case the correlation among modes is due to higher 
derivative operators for single field models with non- 
minimal coupled Lagrangian. Such a correlation is strong 
for modes with comparable wavelength so that the signal 
is maximal for equilateral configurations. Examples of 
models of this kind are the ghost inflation (36| and the 
Dirac-Born-Infcld (DBI) [37f models among others. 

In the following we will focus on the the first class of 
models where non-Gaussianity can be parametrized as 

$(x) = $ £ (x) + / W£ ($|(x) - <$i(x)», (2) 

where $(x) denotes the Bardeen potential, ^(x) de- 
notes the linear Gaussian part of the perturbation and 
the /nl is a merely multiplicative constant that quan- 
tify the level of non-Gaussianity. 

The $(k)-field bispectrum will thus have contributions 
of the form: 



($ L (k 1 )$ L (k 2 )$ JVL (k 3 )) = 
= 2f NL (2 7 r) 3 P$(fc 1 )P (& (fc 2 )<5 3 (k 1 + k 2 + k 3 ) + eye 



(3) 



(i.e. the function F of Eq. [T] is F(fci,fc 2 ,/c 3 ) = 
2fNLP(ki)P(k2) + eye.) where we have used the defi- 
nition of the Bardeen's potential linear power spectrum 4 
P*(A): ($z(ki)$ L (k 2 )> = (2^) 3 P <& (fc 1 ) ( 5 3 (k 1 +k 2 ). 

Note that the non-linearity parameter /jv l defines the 
non-Gaussianity in the gravitational potential and not in 
the CMB temperature fluctuations. 



In order to study higher-order statistics and model 
small deviations from Gaussianity, one can define the 3- 
point correlation function of Bardeen's curvature pertur- 



3 Note that the Newtonian gravitational potential <j> has the oppo- 
site sign of Bardeen's curvature perturbations: <E> = — <f> 

4 Recall that (<t>|(x)) = (2n)~ s J d 3 kP^(k) 



3 



The standard single-field slow-roll inflation model pre- 
dicts a non-Gaussianity of the form described by Eq. @ 

HE mm- 

In standard inflation /tvl is unmesurably small (less 
than 10 -6 [H, Within this picture, the 3-point cor- 
relation function (or equivalently, the bispectrum) turns 
out to be the most sensitive observable to constrain pos- 
sible departures from these (nearly Gaussian) initial con- 
ditions e.g., [2]. Non-linear physics between the end of in- 
flation and the last scattering surface may yield further 
bispectrum contributions (e.g. sec 8 of [pi and references 
therein and more recent work [H, [H, M, EE El ES ED] ) 
which however are expected to be mostly of equilateral 
type and below the detection level for forthcoming CMB 
experiments. A detection of a non-vanishing CMB bis- 
pectrum would be then the smoking gun of a different 
scenario describing the mechanism responsible for the 
generation of the primordial density perturbations. 



B. The Primary CMB bispectrum 

Here we summarize the equations that describe how 
the second-order perturbations in the gravitational po- 
tential translate into perturbations of the CMB temper- 
ature, giving rise to a non-vanishing contribution to the 
CMB bispectrum. 

For adiabatic scalar perturbations the primary contri- 
bution to the CMB coefficients can be written as: 



d 3 k 

(2tt) 3 



$(k) 3T ,(fe)y/ m (n), (4) 



where gTi(k) is the radiation transfer function and 
$(k) is the primordial curvature perturbation in Fourier 
space. From this equation it is clear that, if any, non- 
Gaussianity in <I>(k) will appear in the af m - Accord- 
ing to Eq. ([2]), we can decompose the curvature per- 
turbation into a linear and non-linear term: $(k) = 
$L(k) + $jvi(k) and, by using an analogous notation, 
we will have: a im = af m + af^ . 

Following the steps outlined in Appendix [Al the pri- 
mary CMB angular bispectrum takes the form (l8j : 

/>oo 

C/JT 3(P) = 2 C™ / r 2 dr[bi(r)bl(r)b^(r) 
'Jo 

+ bl(r)bl h (r)bl(r) +&^(r)^(r)^ (r)] 

where Qj i \ m f 2rri3 defines the Gaunt integral (see eq. 
(fATOl) ) and' 3 



k 2 dkP^k)g Tt {k)ji{kr), (5) 



/o 

2 1,00 



b\{r) 

b^ L (r) ee - I k 2 dkh L g Ti (k) Jt (kr), (6) 

JO 

with ji(kr) being the spherical Bessel functions. It is 
important to note that this formula is valid only when 



Jnl does not depend on the scale and it approximately 
applies if such a dependence is weak. Note that for our 
present purposes, if extra contributions are summed to 
the primordial one, we can reinterpret /jyx an effective 
/at l and use the same expression (Eq. [5j for the primary 
bispectrum. Extra contributions are not guaranteed to 
be scale independent or to have exactly the local form, 
making the effective /jvl shape and scale dependent, as 
we will see below. 

Finally, it is useful to define the primary reduced 
bispectrum factorizing the /jyx parameter: &^ 2 ^ = 

Jnl ^i 2 £ 3 j where the quantity frj^ 2 ^ 3 is the reduced bis- 
pectrum for Jnl = 1: 



B 



m 1 m2m 3 (P) 
£1(2(3 



(G mi 



(7) 



C. Secondary bispectra: the cross correlation 
between Lensing and the RS effect 

The path of the CMB photons traveling from the last 
scattering surface can be modifed by the gravitational 
fluctuations along the line-of-sight in several different 
ways. On angular scales much larger than arcminute 
scale the photon's geodesic is deflected by gravitational 
lensing and late time decay of the gravitational poten- 
tial and non-linear growth induce secondary anisotropics 
known respectively as the Integrated Sachs Wolfe (ISW) 
[H[ and the Rees-Sciama (RS) effect [H. Hereafter by 
RS we refer to the combined contribution of linear and 
non-linear growth. 

In this work we will concentrate on the cross corre- 
lation of the CMB lensing signal with the secondary 
anisotropies arising from both the linear ISW and the 
Rees-Sciama effect. We will refer to this as L-RS bispec- 
trum. A closely related effect was investigated in Tlj . 
where only the (linear) ISW contribution is included. Af- 
ter galactic foregrounds, point sources and the Sunyaev- 
Zeldovich [|(| signature from galaxy clusters are expected 
to be the dominant contribution to the CMB bispectrum, 
but because of their frequency dependence and their sta- 
tistical properties they can be separated out without ma- 
jor loss of information [l8|, |4l|. The next leading sec- 
ondary bispectrum contribution is the L-RS one, which 
cannot be separated out by frequency dependence. Both, 
lensing and the RS effect, are in fact related to the grav- 
itational potential and thus are correlated, leading to a 
non-vanishing bispectrum signal with a blackbody spec- 
trum. 

As already pointed out in previous works [H, 0, [HI] 
the joined study of these phenomena through the CMB 
bispectrum is a very powerful tool, for example, to better 
understand linear and non-linear growth of structures, to 
break degeneracies between parameters arising in a power 
spectrum only analysis, or to possibly constrain Dark 
Energy equation of state or models beyond the standard 
ACDM. 
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FIG. 1: TOP-LEFT panel: The non-linear matter power spectrum P^ L (k) obtained with Halofit (solid line) and by using 
the Peacock & Dodds (PD) semi-analytical approach (dot-dashed line). The upper curves refer to redshift z = 0.1, while the 
lower curves to z — 1. BOTTOM-LEFT panel: The absolute value of the Q(£) L-RS bispectrum coefficients defined in Eq.((T3}, 
plotted as a function of the angular scale I. The cusp indicates where Q changes sign due to the onset of non-linearities; in linear 
theory Q is always positive (dashed line). The solid line corresponds to the coefficients obtained by using the Halofit non-linear 
matter power spectra P£ L (k, z), while the dot-dashed line refers to the Q(£) obtained with the PD semi-analytical method to 
model the non-linear behavior. The cosmological parameters used are listed in Tab. [I] Note that the non-linear transition in 
the two cases happens at different scales: the Q(t) from Halofit change sign at I ~ 210, while the ones from the PD at I ~ 300. 
BOTTOM-RIGHT panel: Signal-to-noise ratio -Eq. (|15[l - for the secondary Lensing-Rees Sciama Bispectrum as a function of 
tmax in the case of an all sky, cosmic variance limited experiment (solid line). The dashed line refers to the signal-to-noise for 
the Lensing- linear Integrated Sachs Wolfe bispectrum. TOP-RIGHT panel: The dot-dashed line is the \ 2 between the L-RS 
bispectra obtained respectively with Halofit and with PD, as defined in Eq. 1171 The dashed line represents the same quantity, 
but now the comparison is between the L-RS (Halofit) and the Lensing- linear ISW bispectra. Both quantities are plotted as 
a function of the maximum multipole imax- 



Following [7[ , the CMB anisotropy in a direction h can 
then be decomposed into: 

e{h) = e p {h) + e L {h) + <d RS (h) (8) 

where P denotes primary, L lensing (see Eq. IB1[) and 
RS ISW+Rees-Sciama, which includes both the linear 
and the non-linear contributions. This last term takes 
the form: 

e RS (n) = 2 J dr^4>(r,nr), (9) 

where <j) refers to the gravitational potential perturbation 
and r is the conformal distance defined in Eq. IB31 



We can thus write the bispectrum (see appendix [X| as 

KZT = WCO = ( 10 ) 

(a™ lP a™ 2L a™ 3jR5 ) + 5 Permutations, . 
Following the steps outlined in Appendix IBIthis becomes: 

r,mim 2 m 3 (L — RS) (~>mim2m 3 ,L — RS fTT \ 

where the reduced bispectrum is given by 

,(L-RS) _ ll(ll + l)-l 2 (l 2 + l)+l 3 (l 3 + l) 

x C£ fi(40 + 5 Perm., ( 12 ) 
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and Cf is the primary angular CMB power spectrum. 
Here the quantity that contains physical information 
about the late universe is 0, 0, E3] : 



s r{zi s ) ~ r(z) 
r(z ls )r(z) 3 



—P^ L (k, z) 



(13) 

dz 



that expresses the statistical expectation of the correla- 
tion between the lensing and the RS effect. Here the 
Limber approximation has been used, which we find to 
be extremely good (better than 20%) even at low i's. 
The accuracy of this equation has been explored in [43l |. 
The same coefficients can be calculated in the linear case 
for the cross correlation lensing-Integrated Sachs- Wolfe 
effect (ISW) by simply substituting the non-linear power 
spectrum in equation (|13p with the linear one, P%(k, z). 

In the bottom-left panel of Fig. ([T]) we show the be- 
havior of the absolute value of these coefficients |Q(*)| 
for £ up to 1000 (see 31111 for details). The cusp indi- 
cates that Q(£) changes sign: this is due to the onset of 
non-linearities, which change the sign of dP^/dz. Note 
that in linear theory (dashed line) such a derivative never 
changes sign giving Q(£) > in the A-dominated regime 
(Integrated Sachs Wolfe effect) Therefore this fea- 

ture is a fingerprint of the non-linear regime behavior. 
The scale at which Q(£) changes sign depends crucially 
on the scale at which the non-linear growth overcomes 
the linear effect, making the the L-RS bispectrum sensi- 
tive to cosmological paramete rs g overning the growth of 
structure like Q m , w or er 8 [71. Ha. |20I| . 



III. BISPECTRUM CALCULATION AND 
EXPECTED SIGNAL-TO-NOISE 

We assume a fiducial ACDM model in agreement with 
the latest observational results [24j with the parameters 
listed in Tab. HI The L-RS bispectrum calculation, Eq. 
(TT21 [T3"l) . requires evaluation of the non linear P ( j > (k,z). 
The gravitational potential 4> is related to the matter 
density fluctuation 6 through the Poisson equation: 

P*(M) =("«.») (jf) P S (k,z)(l + z) 2 . (14) 

where f2 m = fif, + f2 c is the total matter density parame- 
ter. There are two approaches to compute the Pg fL (k, z) 
that have been extensively tested and used in the liter- 
ature: the more recent Halofit [44( model, which is in- 
cluded in CAMB [H[ and the Peacock and Dodds (PD) 
46 1 method (generalized for dark energy cosmologies by 
47|). Here we use both approaches and compare them. 
Note that the literature so far on the L-RS bispectrum 
0, 0, H3| has used the PD approach to describe non- 
linearities. 

We perform numerical derivatives to map the function 
dP£ rL (k,z)/dz at k = -JU with £ up to 1500. Then, to 



TABLE I: ACDM parameters. 



Symbol 


Description 


Value 




Hubble constant 


70 Km/sec/Mpc 


n h 


Baryon density 


0.044 




Dark matter density 


0.224 




Dark energy density 


0.732 


w 


Dark energy equation of state 


-1 




Fluctuation amplitude at 8/i _1 Mpc 


0.834 


n s 


Scalar spectral index 


1 


Zls 


Redshift of decoupling 


1090.51 



compute the Q(£) coefficients (see bottom-left panel of 
Fig. (HJ)) , wc numerically integrate in < z < 2.5, this 
is sufficient to account for the dark energy signature and 
the non-linear regime (widening the integration interval 
does not change the results) . For the primary bispec- 
trum, we proceed as in [g, [18| assuming a ACDM model. 
We compute the radiations transfer functions gTt(k) with 
the CMBFAST code [H[ and weperform the k and re- 
integrations in the same way [1, [l8| did. 



A. The Lensing-Rees-Sciama bispectrum: 
Signal-to-Noise ratio 

According to 4], the bispectrum signal-to-noise ratio 
can be, in general, defined as: 



iii2i3 



(Bp 



(15) 



where A^,,^ is a number which takes value 6 for equi- 
lateral configurations, 2 for isosceles configurations and 
1 otherwise (see [49( for details). 

Using Eq. (|A12[) we can write the the signal in the 
numerator as: 



(Bo 



(24 + 1) (2& + 1) (2£ 3 + l) 



4tt 



(16) 



OOO; "'I***' 



In the bottom-right panel of Fig. Q]the solid line shows 
the signal-to-noise ratio for the L-RS bispectrum as a 
function of the maximum multipole £ max [£\ , £2 and £3 
are all < £ m ax)- The dashed line refers to the signal-to- 
noise for the linear case (L-ISW bispectrum). Note the 
enhancement due to non-linearities at high multipoles. 
We do not consider £ ma x > 1500 because other secondary 
effects (e.g., Ostriker-Vishniac or Kinetic SZ jSHSHJ) 
may start to dominate. The S/N plotted has been ob- 
tained by summing over all triangle configurations for 
a full sky, ideal, cosmic variance-dominated experiment. 
The results can be representative of an experiment with 
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the nominal performance of Planck, as pointed out in 
previous works [2(| H2] . 

The signal-to-noise ratio increases mainly when the 
maximum multipole £ max reaches few hundred, where 
the signal gives the main contribution. As we will ex- 
plore in more detail later on, the L-RS bispectrum signal 
dominates for squeezed triangle configurations when a 
large scale mode couples with two small scales modes: 
50% of the signal-to-noise comes from triangles with 
2 < (-mm < 10, in agreement with the findings of [tJ- 



B. Modeling non-linearities: Peacock & Dodds and 
the Halofit model 



The two main approaches that can be used to compute 
the non-linear matter power spectrum Pg fL (k) are the 
commonly used semi-analytical Peacock & Dodds [46[ 
(PD) formula, based on the scaling method of [H[ , and 
the more recent Halofit model |44j |. 

The first approach is based on the ansatz that the non- 
linear evolution induce a change of scale so that the non- 
linear power spectrum at wavenumber k can be param- 
eterized by a simple function of the linear one evaluated 
at k 1 . This has been shown to interpolate correctly the 
P(k) behavior in the intermediate regime between linear 
and stable clustering. 

The second approach is based on the so-called "halo 
model" for the matter power spectrum. In the halo 
model the density field is decomposed into a distribution 
of clumps of matter with some density profile. The large 
scale behavior is then derived through the correlations 
between different haloes, while the non-linear correlation 
functions on small scales are obtained from the convolu- 
tion of the density profile of the halo with itself. Halofit 
has been also extensively tested on large, high-resolution 
N-body simulations. 

We estimate that any uncertainty in the description of 
the non-linear clustering should be at or below the level 
of the difference between these two approaches. 

In the top-left panel of Fig. [1] the non-linear matter 
power spectrum Pg !L (k) is plotted as a function of the 
wavenumber k for Halofit (solid line) and for PD (dot- 
dashed line) . The upper curves refer to power spectra at 
redshift z — 0.1, while the lower curves are the non- linear 
matter power spectra at z = 1. The Halofit power spec- 
trum shows the Baryon Acustic Oscillation (BAO) at the 
typical BAO scale k ~ O.lMpc . To produce the PD one 
we started from a "no-wiggle" linear power spectrum. 
This is because the PD approach maps linear scales into 
non-linear ones thus artificially changes the position of 
the wiggles; when taking derivatives this can induce spu- 
rious signal which does not happen when staring from a 
"no- wiggle" linear P(k). Beside the BAO feature, which 
is irrelevant for our purpose, the two models are in good 
agreement although at higher z the PD power spectrum 
seems to produce a power spectrum more non-linear that 



Halofit. 

The bottom-left panel of Fig. [1] show the effect of 
this difference in the L-RS Bispectrum coefficients Q(£). 
The figure shows the absolute value of the coefficients 
Q(^)| as a function of the angular scale I. The solid line 
corresponds to the coefficients obtained by using Halofit 
while the dot-dashed line using PD. The transition to 
the non-linear regime (indicated by the cusp where Q(£) 
changes sign) happens at smaller £ for Halofit (~ 200) 
than for PD case (I ~ 300). 

We quantify the difference between the two models by 
computing the % 2 for the L-RS bispectra obtained re- 
spectively with Halofit and with PD: 



XHalofit-PhD 



E 



( „L-RS [Halofit] _ „L-RS[P&cD]\'' 



(17) 

This is shown in the top-right panel of Fig. Q] (dot-dashed 
line) where it is plotted as a function of the maximum 
multipole £ max for our fiducial cosmology. The two mod- 
els are compatible within 1-cr (A^ 2 < 1) for £ max < 900. 
We conclude that the significance of a detection of this 
signal does not depend crucially on the modeling for non- 
linearities, however the choice of an incorrect modeling 
may introduce significant biases when doing precise anal- 
ysis, as 0,13) on key parameters, e.g., w, as, fi m which 
are particularly sensitive to the onset of non-linearity. 

For each of the parameters w, as, Q m we compute the 
bias introduced by using Halofit in the bispectrum cal- 
culation in the hypothetical case that PD was a true de- 
scription of non-linearities. Around our fiducial model, 
the biases are at the level comparable to the 1 — a er- 
rors (0.3 to 0.6 a). We estimate that, in any practical 
application, biases introduced by uncertainties in the de- 
scription of non-linear clustering will be at this level or 
below. Moreover, dP NL /dz can be accurately evaluated 
with the use of N-body simulations as presented in [54| , 
thus removing this source of bias. Ultimately, in the top- 
right panel of Fig. Q]we plotted the x 2 from the L-RS and 
the linear Lensing-ISW bispectra. This is defined as in 
Eq. (fT7|) . but with the Lensing-ISW bispectrum instead 
of the non- linear one obtained by using the PD model. At 
high multipoles the linear and the non-linear cases differ 
by more than 1-a. The error introduced by the change 
in modeling non-linearities is smaller with respect to the 
error due to only consideing the linear behavior. 



IV. THE SHAPE-DEPENDENCE OF THE 
BISPECTRUM SIGNALS: PRIMARY VS L-RS 

The signal for the primary bispectrum is dominated 
by squeezed and nearly squeezed configurations as shown 
by [H, [H, H^] . We find that the same applies to the L- 
RS bispectrum, where 90% of the signal-to-noise comes 
form nearly squeezed configurations where £2 > 10£i and 
£■2 < £3- This can lead to contamination, i.e. "confu- 
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200 400 600 800 1000 200 400 600 800 1000 

l z >40 (L 3 =l 2 +2, L,=2) l z >40 (L 3 =l s +2, L,=2) 




FIG. 2: Effective non-linear parameter /j^£ flS (Eq. I18I left panels) and corresponding reduced L-RS (solid) and Primary 
(Dashed) bispectra (right panels) for two nearly squeezed configurations: £\ = 2, £2 > 40, £3 = £2 + 2 (top panels) and 
£\ = 2, £2 > 40, £3 = £2 (bottom panels). In the right panels the primary bispectrum plotted has /jvl = — 10 for making it 
more visible. Note that in the case illustrated in the bottom panels, the two bispectra have exactly the same shape and they 
are completely degenerate for a /nl — —17. 



sion" , between the two signals. We illustrate this point 
by defining an "effective" Jnl for the L-RS signal as 

iL-RS 

UNL It-Ms - JT, > (J-oj 

which depends on triangle shape. This is shown in the 
left panels of Fig. [5] while the right panels show the cor- 
responding reduced bispectra (solid for L-RS and dashed 
for primary) . The top panels are for nearly-squeezed con- 
figurations where £\ is fixed, l\ = 2, £ 2 varies for £ 2 > 40, 
and £3 = £2 + 2 while the bottom panels are for isosceles 
squeezed configurations: £\ = 2, £ 2 > 40, £3 = £2- Note 
that in the right panels the primary bispectrum has been 
computed for /ml = — 10 for making it more visible. 

The case of the squeezed isosceles configurations 
(where £2 ^> £1, £3 — £2 and £\ <~ 150), contributes 
with a ~ 5% to the total S/N in both cases, and the 
two bispectra have exactly the same shape and they 
completely degenerate for a f^L — — 17- For nearly 
squeezed-configurations f^1 RS oscillates but its average 



is at around /jvl^ 5 ~ 10. Nearly squeezed configura- 
tions where £ 2 > 10£i and £ 2 < £3 carry most of the S/N. 
The L-RS signal always dominates over the primary one, 
normalized to /jvl = 1, for all the high signal-to-noise 
configurations by a factor of 10-20 in absolute value. 

However, besides to the fact that the two bispectra 
could be confused for showing some similar behavior (see 
bottom panels of Fig. [2]) they can in principle be dis- 
entangled since they have intrinsically different features 
arising from the extremely different physics behind them. 

For example, looking back at the top-right panel of 
Fig. [21 we find that for these configurations the L-RS sig- 
nal oscillate, while the primary reduced bispectrum does 
not. The L-RS bispectrum of eq. (IT2|) in fact contains 
the Ci, with the typical structure given by the acoustic 
peaks, and the coefficient Qp which determine the change 
of sign. On the other hand, the primary signal, see eq. 
(O, is composed by the coefficients: bf(r) oc P^,(k)gTe(k) 
and bf L (r) oc fNLgre(k) so that the changing of sign in 
this case is due to the full radiation transfer functions 
gxeik). For general configurations, the two bispectra be- 
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FIG. 3: Reduced bispectra bi ± i 2 i 3 : Primary for /„( = —10 (dashed line) and Lensing-Rees Sciama (solid line). The left plot 
shows equilateral triangle configurations I = i\ = 1% = 1%, while the right one shows the flattened configurations £\ — 2£s and 
li = £3- We have plotted bf_J^ £\(£\ + l) 2 10 16 , which makes the Sachs-Wolfe plateau of the Primary reduced bispectrum 
easily seen at large angular scales. 



have differently: in Fig. [3] we plot the case of equilat- 
eral (left panel) and flattened configurations of the type: 
l\ = 2^3 and £2 = I3 (right panel). The dashed lines 
refers to the primary contribution while the solid lines to 
the L-RS one. The two bispectra have different shapes 
and change sign at different angular scales. In the case of 
equilateral configurations, for example, the primary re- 
duced bis pec trum shows the known oscillatory shape, as 
found in [Tg|, while the L-RS reduced bispectrum does 
not. In the case of flattened configurations both bispec- 
tra show oscillations. Note that in these plots the y-axis 
has been multiplied by a factor 10 16 (while in Fig. |(3J) the 
y-axis has been multiplied by a factor 10 11 ): these con- 
tributions are clearly subdominant by about 5 orders of 
magnitude with respect to the squeezed configurations, 
which explains why the latter shapes dominate the signal- 
to-noise. 

In light of these findings we now attempt to interpret 
recent constraints on prim ordial non-Gaussianity from 
CMB data e.g., [2TI. [23ll2^| and consider the implications 
for forthcoming measurements. 

The $nl estimator used in these works reduces to the 
one defined in [56] in the simplest case of temperature- 
only anisotropics, cosmic variance dominated, all sky 
analysis. This estmator weights the bispectrum of every 
triplet ^i,^2)^3 by the signal to noise of the primary bis- 
pectrum. We can thus estimate the contamination that 
such an estimator would measure due to the presence of 
the L-RS signal defining: 



Jnl 



S_ 

N 1 



where 



and 



E 



D 



B 



L-RS 

Ce 1 C( 2 Ci 3 



N 



E 



Cf_, C( 2 Ci 3 



(20) 



(21) 



(19) 



This is plotted as a function of l ma x in Fig. [51 up 
to £ ma x = 1500: the dashed line refers to Jnl ob- 
tained by summing over all configurations, while the 
dot-dashed line refers to Jnl obtained from only nearly 
squeezed configurations (X^°= 2 Y.^oii T,e7=W' which 
dominate for both the primary (local type) and the 
lcnsing-Rees Sciama bispectrum. The solid lines indi- 
cates where the bias is negative. This is in qualitative 
agreement with the effect explored in illj where only lin- 
ear growth was included: the inclusion of (later-type) 
non-linearities enhances the contamination but, as ex- 
pected, only for I > 400. For example at £ ~ 1500 we 
estimate an enhancement of roughly 40%. 

Using Halofit or PD in the modeling of the L-RS bis- 
pectrum changes the estimates of the effective Jnl by 
10% indicating that this correction is robust to possible 
residual uncertainty in the modeling of non-linearities. 

It is also possible to estimate the effective /nl via a 
simple x 2 analysis: we find the same values as above but 
the interpretation of the \ 2 as a goodness of fit test would 
indicate that the local model is not a good fit. This is in 
qualitative agreement with the findings of [57| who com- 
pute the CMB bispectrum from the second-order fluctu- 
ations and find that their effect is separable from the pri- 
mary non-Gaussian signal because of the different shape 
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dependence for non-squeezed (or nearly squeezed) con- 
figurations. The agreement cannot be made fully quanti- 
tative as perturbation theory approach may break down: 
for a given multipole £ the derivative of the gravitational 
potential power spectrum is probed at a wide range of 
scales k(z) — and therefore highly non-linear scales 
can contribute non-negligibly even at relatively low I. In 
practice, however, it may not always be possible to im- 
plement a goodness of fit test. 

The expected error on /jvl for forthcoming surveys is 
smaller than 10 (for example the Planck surveyor, re- 
cently launched is expected to yield 1-a error on /jyi, 
of order 4 52]), indicating that the L-RS signal may be 
a crucial contaminant in the pursuit of primordial non- 
Gaussianity, if not properly taken into account. We have 
shown here that its amplitude and configuration depen- 
dence is well known; it is thus not necessary to extract 
this signal from the CMB bispectrum and separate it 
from the primary: it can simply be included in the mod- 
eling of the CMB bispectrum. 



V. CONCLUSIONS 

We have revisited the predictions for the expected 
CMB bispectrum signature of the primary-lensing Rees- 
Sciama (L-RS) correlation. This bispectrum is the lead- 
ing secondary contribution on scales much larger than 
arcminute with the same frequency dependence as the 
CMB primary. Forthcoming experiments like Planck 
have the statistical power to detect this signal with 
a signal-to-noise of order 10. The linear contribution 
(primary- lensing-ISW bispectrum) was considered in 
H, [|[ and in [ll[ ■ By including the non- linear (RS) de- 
scription the signal-to-noise increases to ~ 10 to £ m ax = 
1000. The overall signal depends on the balance of two 
competing contributions along the line of sight: the de- 
caying gravitational potential fluctuations and the ampli- 
fication due to non-linear gravity. For this reason the ef- 
fect can be used to place strong constraints on cosmolog- 
ical parameters that determine the growth of structures: 
f2 m , dark energy parameters and a%. By comparing two 
different semi-analitic descriptions of non-linear cluster- 
ing, we find that an accurate description of the non-linear 
growth of the matter power spectrum is necessary to ob- 
tain unbiased estimates of these parameters. Approaches 
based on numerical simulations (see e.g., [HI, [Hj]) will 
have to be employed. In general, while the approxima- 
tions used here to derive and compute Eq |13l are ex- 
tremely good for the purpose of this paper, a detailed 
comparison with data will require the exact numerical 
evaluations. 

Here we have shown that this bispectrum signal can 
be confused with the signal from local primordial non- 
Gaussianity. Both bispectra signal are maximal for 
squeezed or nearly squeezed configurations. For some 
configurations (e.g., squeezed isosceles) the two bispec- 
tra are virtually identical, while for generic configurations 
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FIG. 4: The plot shows }nl, as denned in Eq. I19I as a func- 
tion of Imax- The dashed line refers to Jnl obtained by sum- 
ming over all configurations, while the dot-dashed line refers 
to /ml obtained from nearly squeezed configuration (where 
£i runs from 2 to 10, £i from 50^i to £ ma x and £3 from 1% to 
£max), which dominate for both the primary (local type) and 
the lensing-Rees Sciama bispectrum. The solid lines indicates 
where the bias is negative. 



the shape dependence of the two bispectra are different in 
the details. A bispectrum estimator optimized for con- 
straining primordial non-Gaussianity of the local type 
would measure an effective Jnl — 10 for £ max — 1000 due 
to the presence of the primary- lensing-Rees-Sciama cor- 
relation. If not accounted for, this introduces a contami- 
nation in the constraints on primordial non-Gaussianity 
from the CMB bispectrum. This is in qualitative agree- 
ment with the effect explored in [ll| where only linear 
growth was included. For £ > 400 the full non-linear 
treatment is needed. For current data, this contamina- 
tion (effectively bias in the recovered /jvx) is smaller than 
the 1-a error, however it can become significant when in- 
terpreting the statistical significance of results that are at 
the boundary of the 3 — a confidence level. For example 
if we subtract the effective value for the L-RS Jnl from 
the central value of the estimate of [2l[ , we obtain that 
/nl primordial is consistent to zero at the ~ 2.5a confi- 
dence level. For forthcoming data, however, this bias will 
be larger than the 1 — a error and thus non-negligible. 
A more quantitative statement cannot be made at this 
stage because the calculations presented here are done 
for a cosmic-variance dominated experiments while for 
the current bispectrum analysis from WMAP data in- 
strumental noise cannot be neglected at £ ~ 400 where 
most of the contamination is expected to come from. 

Techniques to separate out different bispectra shapes 
and assess whether a detection of non-Gaussianity is pri- 
mordial have been proposed [HI] and will be suitable for 
this application. 

We argue that the bispectrum of the L-RS effect can 
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be accurately modeled: even with currently available 
semi-analytic descriptions for non-linear clustering, we 
estimate the error on the effective Jnl to be at the 
10% level or below. We conclude that, in analyzing 
the CMB bispectrum to obtain constraints on primordial 
non-Gaussianity for forthcoming data, this contribution 
must be included in the modeling. 
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APPENDIX A: BISPECTRUM STATISTICS 

Deviations form Gaussianity in the CMB are charac- 
terized by the angular n-points correlation function of 
the temperature field in the sky 



(e(n 1 )6(n 2 )...e(n n )) 



(Al) 



where the bracket defines the ensemble average and n 
the angular position (i.e. the direction unit vector of the 
incoming photons). In general it is useful to expand the 
field in terms of spherical harmonics: 



8(n) = ^2 a im Y hn (h), 



(A2) 



e=o m=- 



so that, by using the symmetric proprieties of harmonic 
transformations and the orthogonality of the spherical 
harmonics, we can write the coefficients a" 1 as: 

af = J d 2 n e(n)r/ m (n). (A3) 

The angular CMB bispectrum is defined by three har- 
monic transforms satisfying rotational invariance: 

thus the angular averaged bispectrum takes the form: 



B 



E 



all 



4 4 4 



ii m 2 m 3 J Mats 



B): 



(A5) 



Since 4 5 4 and 4 form a triangle, this quantity must 
satisfy the triangle conditions and parity invariance: 

mi + m 2 + m 3 = 0, l\ + 4 + 4 — even, 



(A6) 



for all permutations of indices. The matrix appearing in 
equation (jA5|) represents the Wigner-3j symbol that de- 
scribes the coupling of two angular momenta. Rotational 
invariance requires the bispectrum amplitude to be inde- 
pendent from orientation and triangle configuration. The 
Wigner-3j symbol, transforming the m's under rotations, 
preserve the triangle configuration thus describing the 
bispectrum azimuthal angle dependence. The orthogo- 
nality properties of the Wigner-3j symbols are 



E 



E 

all in 



4 4 4 

m'i m'n m'-i 



4 4 4 

mi TO2 TO3 



= 1 



(A7) 



4 4 L 
m[ m' 2 M' 



'M' 



2L + 1 



(A8) 

By making use again of rotational invariance and of the 
simmetry and orto-normality properties of the 3-j sym- 
bols, we can write the bispectrum as: 



r>m 1 m 2 m; i 



r >m 1 m 2 in 3 i 
Villus "iiWs 



(A9) 



where G^l^" 13 is the Gaunt integral which contains all 
the angle dependence and triangle constraint information 
and it is defined by: 



J d 2 nY eimi (h)Y t2m2 (n)l^ 3m3 (n) 



(2£i + 1) (2£ 2 + 1) (2*3 + 1) 



4tt 



(A10) 



4 4 4 





4 4 4 

mi m 2 7713 



where b^i^ is called the reduced bispectrum, which is 
a very useful quantity since it is an arbitrary symmet- 
ric function of 4,4 and £3 only and it contains all the 
relevant physical information of the bispectrum. 

By substituting equation (|A9[) into equation (|A5|) and 
using the Gaunt integral property: 



E 

all m 



t\ £2 £3 \ rtmim^m-^ 

mi m<2 7713 / £ i £ 2^3 



(All) 



(24 + 1) (24 + 1) (24 + 1) fli £ 2 £3 





4tt 



we can finally write: 



B 



(24 + 1) (24 + 1) (24 + 1) 



-in 



(A12) 



4 4 4 
000 
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For high-f the Gosper factorials approximation for the 
Wigner 3j symbols can be used: 



h h h 




L 



L+l 



L/2 



1 



'■->< 1L- 'V /2 -q (6L- 12^ + 1) 



7T L+l 



(6L + 7) 1/4 

1/4 



T7i( A13 ) 



1/2 



APPENDIX B: WEAK LENSING OF THE CMB 

Weak lensing of the CMB re-maps the temperature 
primary anisotropy according to: 



6 L (n) = 9 p (n + V(/>) 

~ 6 p (n) + Vi<Kn)V'e p (n) 



(Bl) 



where the label 'L' refers to the lensed term while 'P' 
to the primary contribution. The deflection angle a = 
V</>l is given by the angular gradient of the gravitational 
potential projection along the line of sight: 

cb L (h) =-2 [ ' 3 dr r( - ZLs) ~ T{Z \ (r,nr). (B2) 
Jo r{z)r{zis) 

Here r is the comoving conformal distance. Assuming 



a flat ACDM universe this can be written as: 

dz' 



r(z) 



(B3) 



and thus ri s = r(zi s ) refers to the comoving radius at 
last scattering from the observer at z = 0. 

As done with the temperature perturbations, we can 
expand the lensing potential into multipole moments: 



(B4) 



By applying eq. (|A3[) into eq. (|Bip and carrying out 
the calculations we get an explicit expression for the lens- 
ing a™ coefhcients: 



aT L = aT P 



E c-i: 



m+m +m r>—mm m /-pc\ 



£'(£' + !)-£{£ +!)+£" {£" + !) m , p , „ 
x - a e , (p Le „ 



being C?^,™V m the Gaunt integral defined in equa- 
tion(fST0|. 
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